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The primordial non-Gaussian parameter /nl has been shown to be scale-dependent in several models of 
inflation with a variable speed of sound, such as Dirac-Born-Infeld (DBI) models. Starting from a simple ansatz 
0^ for a scale-dependent amplitude of the primordial curvature bispectrum for two common phenomenological 

models of primordial non-Gaussianity, we perform a Fisher matrix analysis of the bispectra of the temperature 
and polarization of the Cosmic Microwave Background (CMB) radiation and derive the expected constraints on 
the parameter ;jng that quantifies the running of /nl(^) for current and future CMB missions such as WMAP, 
(— I Planck and CMBPol. We find that CMB information alone, in the event of a significant detection of the non- 

r-j Gaussian component, coiTesponding to /nl = 50 for the local model and /ml =100 for the equilateral model of 

K— 5 non-Gaussianity, is able to determine «ng with a 1-cr uncertainty of A«ng - 0.1 and Awng - 0.3, respectively, 

I for the Planck mission and a factor of two better for CMBPol. In addition, we consider a simple Fisher matrix 

analysis of the galaxy power spectrum to determine the expected constraints on the running parameter wng 
for the local model and of the galaxy bispectrum for the equilateral model from future planned and proposed 
photometric and spectroscopic surveys. We find that, in both cases, large-scale structure observations should 
achieve results comparable to or even better than those from the CMB, while showing some complementarity 
due to the different distribution of the non-Gaussian signal over the relevant range of scales. Finally, we compare 
our findings to the predictions on the amplitude and running of non-Gaussianity of DBI inflation, showing how 
^ I the constraints on a scale-dependent /nl(^) translate into constraints on the parameter space of the theory. 
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I. INTRODUCTION 

Present observations still have little power in discriminating between the huge number of different inflationary sce- 
narios that have been proposed so far in the literature. Moreover, alternatives to inflation such as cyclic and ekpyrotic 
^y-^ Universes are also compatible with the data. In order to discriminate between all these possibilities we need to move 

from the study of general properties shared by all the models under examination, such as flatness and scale-invariance, 
C3 to more specific and strongly model-dependent predictions. For example, diff'erent inflationary models predict the 

\^ primordial curvature perturbations to be close to but not exactly Gaussian. The specific departure from Gaussianity is 

highly model-dependent and, if detected, can be an extremely valuable tool for the purpose of discriminating between 
alternative scenarios. 

Until fairly recently it has been generally believed that measurements of the bispectrum of CMB anisotropics, |[T]|2] 
jj^ IS 4, 5 1, were able to put the most stringent constraints on /nl when compared to all the other alternatives. However, 

the use of large-scale structure (LSS) observations as a probe of primordial non-Gaussianity has received increasing 
attention in the last year as it has been shown that present and future galaxy surveys allow one to obtain constraints 
at the same level of those from the CMB. Non-Gaussian initial conditions are expected to have a significant effect on 
the high-mass tail of the halo distribution, (|6]|7l|8l|9l[T0l[lTl[T2][l3], and it has been more recently recognized in the 
void distribution Iil4l . A large amplitude of the primordial bispectrum could lead to an observable initial component 
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in the large-scale skewness of the galaxy distribution, ifTSl [161 [TtII or on the galaxy bispectrum, ifTSl [191 [20l |2T]. 
More recently, significant interest has been generated by the unexpected but large effect of local non-Gaussian initial 
conditions on the bias of halos and galaxies |[22l [23l l24l 25 26 27] |28] (29 1. Remarkably, this new method lead 
to constraints on the specific local type of non-Gaussianity from current data-sets, which are already competitive 
with CMB limits, |l24','26|. Finally, additional confirmation of a possible detection might also come from alternative 
statistics such as Minkowski Functionals ll30l[3Tl . 

So far, measurements of the bispectrum of the CMB temperature fluctuations, a direct probe of the primordial 
bispectrum, have been found to be consistent with the Gaussian hypothesis and limits have been placed on the am- 
plitude of a possible non-Gaussian component, ||5]|32l|33]. Such limits are placed on the amplitude of the curvature 
bispectrum — the /nl parameter — assuming the specific dependence on the shape of the triangular configuration pre- 
dicted by the inflationary model. This amplitude is usually assumed to be scale-independent. However, models such 
as Dirac-Bom-Infeld (DBI) inflation |[34ll35l[36l . can be characterized by a primordial bispectrum whose amplitude 
varies significantly over the range of scales accessible by current cosmological probes. 

Motivated by these considerations, the aim of this paper is to study what limits can be placed, in the event of a 
detection, on the scale-dependence of the bispectrum amplitude using present and forthcoming cosmological data- 
sets and how such information can be used to obtain additional constraints on inflationary models. This issue was 
previously addressed by Lo Verde et al. |[T3l , who made forecasts based on cluster number counts as a probe of non- 
Gaussianity on small scales combined with current CMB constraints. This work, however, simply assumed previous 
CMB limits on /nl as constraints on the amplitude alone at a given pivot scale. In this regard, we make a step 
forward by providing a complete Fisher matrix analysis of the CMB bispectrum described by two parameters: the 
overall amplitude, /nl, and a non-Gaussian running parameter hng. In addition we combine these results with simple 
estimates of those achievable by LSS observations such as the galaxy power spectrum and bispectrum in forthcoming 
surveys. 

The paper is organized as follows: In Section[ll]we discuss various parameterizations of the primordial bispectrum 
and the phenomenological ansatz assumed in our analysis. In Section III we review the CMB bispectrum, define our 
notation and derive the CMB Fisher matrix for the amplitude /nl and running hng, showing the results for current and 
future CMB missions. In Section IV we derive the Fisher matrices for the LSS power spectrum and bispectrum and 
show the corresponding constraints for a sample of future photometric and spectroscopic galaxy surveys. In SectionjV] 
we present the combined CMB plus LSS results. In Section VI we consider DBI as an example of inflationary model. 
Finally in Section VII we summarize our findings and give some concluding remarks. 



II. PARAMETERIZATIONS OF PRIMORDIAL NON-GAUSSIANITY 



The statistical properties of a Gaussian random field are encoded in its two-point correlation function or, equiva- 
lently, its power spectrum since all higher-order connected correlation functions are vanishing. On the other hand, 
a non-Gaussian field is given, in principle, by an infinite set of functions. Any model of inflation indeed provides 
predictions for all primordial correlators which are characterized not only by an overall constant amplitude but also 
by their peculiar dependence on the shape of the configuration of points in position space or wavenumbers in Fourier 
space. 

In this Section we introduce two of the most common functional forms for the curvature primordial bispectrum, 
the so-called local and equilateral shapes, and describe the specific ansatz for the scale-dependence of the bispectrum 
amplitude that we assume in our analysis. 



A. Shape-dependence and amplitude of the primordial bispectrum 

The non-Gaussian initial conditions can be generally characterized in terms of the bispectrum B(s,{k\,k2,k-i) of 
Bardeen's curvature perturbations defined as 

((Dk,<l>k,<l>k3) = 4(ki23) B^ikukM, (II. 1) 
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where we introduce the notation k,j = k, + kj. We can separate the overall amplitude /nl from the functional form 

F{k],k2,kj,) as 

B^{h,k2,h) = /nl F{hM,h), (11.2) 

where we will define the normalization of F below. Here /nl is a dimensionless parameter while the function 
F(ki,k2,ki) describes the dependence on the shape of triangular configuration defined by the three wavenumbers 
ki, k2 and k^, typically having the hierarchical behavior F{k,k,k) ~ P^(k) for equilateral configurations, with P,s,(k) 
being the curvature power spectrum. In principle, the curvature trispectrum is also expected to have a significant effect 
on large-scale structure observables at the largest scales |37', '38l|. However, for simplicity we will ignore effects due 
to correlations functions of order higher than the bispectrum. A more detailed discussion on their role will be given in 
Section HV] 

Different inflationary models predict different values for /nl, starting from 0{.Q1) up to very large values. CMB 
observations provide the upper limit /sjl < 100, already constraining some of the existing models. A theoretical lower 
bound on the detection of primordial non-Gaussianity is roughly /nl <: 0(1), below which second-order perturbations 
from post-inflationary evolution become relevant. Non-Gaussianity from canonical single-field slow-roll inflation 
models is predicted to be very small O(.Ol) ll39l l40l |4TI l42l l43l l44ll : however, a large class of more general models, 
e.g. models with multiple scalar fields, features in inflation potential, non-adiabatic fluctuations, non-canonical kinetic 
terms, deviations from the Bunch-Davies vacuum, among others, predict substantially higher level of primordial non- 
Gaussianity (for a review, see |45 1 and references therein, for recent contributions see also B6l l47l l48l ). For this 
reason alone, a detection of primordial non-Gaussianity would have the important consequence of ruling out canonical 
single-field slow-roll inflation as a viable inflationary scenario. 

In addition to the amplitude differences, different models of inflation also lead to different functional forms for the 
bispectrum, characterized by F(ki,k2,k-i). Such functional forms can be broadly classified into three classes |,49ll50l 
ISTl : a local, "squeezed," shape-dependence where F{k\,k2,ki) is large for the configurations in which k\ «i k2, k^ 
(and permutations); a non-local, "equilateral," where F{k\,k2,k-i) is large for k\ ~ k2 ~ k^ and a "folded" shape- 
dependence where F{ki,k2,kj,) is large for flattened configurations ki ~ k2 ~ kj,/2 (and permutations). Several 
other shapes, however, have been considered in the literature. For a more complete discussion on shapes and their 
correlations, see for instance Fergusson and Shellard |50|. We limit our attention to the local and equilateral forms 
both because of their simplicity and because they mimic a large part of the models present in the literature. 

The local form can arise from a non-linear relation between the inflaton and curvature perturbations |39, 40, 42], or 
alternatively in curvaton models |52 1. Alternatives to inflation like New Ekpyrotic and cyclic models are also expected 
to produce a large level of non-Gaussianity of this type Il53ll54ll55ll56l . The local form of non-Gaussianity is so-called 
because it can be parametrized in real space by the local expression |fTl [T9ll42l : 

0(x) = (Dg(x) + f^l- [0.2 (X) - <<d2 (x))l , (II.3) 

where is the Gaussian part of the perturbations. The term {0?,(x)) ensures that the perturbation has zero mean. 



From Eq. (II. 3 1 it is easy to show that the shape dependence, F\ac,{ki,k2,k^), for the local model takes the following 
form: 



Fioc.(ki,k2,ki) = 2A| 



1 1 1 



,3-(«,,-l),3-(«,,-l) ,3-(«,,-l),3-(«,-l) ,3-(«,-l),3-(n,-l) 
ft, J ft^ X 3 2 3 



(11.4) 



where the normalized power spectrum is defined in terms of the tilt n, and the curvature power spectrum P<s,(k) 
as P,s,(k) = At^fk^^^^"'^^^ with P(t,{k) defined in terms of the Gaussian component alone as ((I)G(ki)<l>G(k2)) = 

6f(ki2)PM- 

Equilateral forms of non-Gaussianity arise from models with non-canonical kinetic terms such as the DBI ac- 
tion |35J, ghost condensation [57J, or any other single-field models in which the scalar field acquires a low speed of 
sound 1581 \59l . Although the shapes predicted by different models are in this case not identical, it has been noted 
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that they are all very well-approximated by the function: 



r— — — --— — - + 2 perm. 



\l-(«.-l)/3,2-2(n,-l)/3,3-(«,-l) ^ P^™' 



(11.5) 



The defirution for the equilateral model follows from the local one since is defined in such a way that for equilateral 



configurations, F^^ik,k,k) 



F\oc.{k, k, k) and one obtains the same value for Z?o given - f^^. 



Comparisons with observations of different forms of non-Gaussianity lead to constraints on the different amplitude 
parameters. The most recent analysis of the CMB bispectrum provides for the local non-Gaussian parameter the limits 
-4 < < 80 at 95% C.L. corresponding to a 1 cr error on of A/jJjjJJ- ^ 21 [5|. In the equilateral case we have 
-151 < < 253 with A/j^^ ^ 101 ||32|| . Future missions are expected to provide uncertainties of the order of 
~ 3 and Af^l ~ 10, HU 13111. Recent constraints on local non-Gaussianity from the bias of high-redshift objects 
correspond to A/^°j^- ~ 24 ll24l . Future large-scale surveys might probe f^^- with A/^™ ~ 1, ll24l 12911611 . This method, 
however, would lead to very mild constraints for the equilateral form of the initial bispectrum which is expected to 
have a negligible effect on halo bias, ||23ll27l . On the other hand, expected constraints on equilateral non-Gaussianity 
from measurements of the galaxy bispectrum in spectroscopic surveys are about A f^ ~ 25 |21 1. 

Before concluding this Section, we note here that two distinct definitions of /nl are present in the literature, cor- 
responding to a CMB convention and a LSS convention. In the CMB convention, for local non-Gaussianity, /'°'^- is 



defined by Eq. (II. 3 i with the curvature perturbations O evaluated at early times during matter domination, when their 
value its constant. In the LSS convention, one usually assumes O to be the value linearly extrapolated at present time, 
and therefore includes the late-time effect of the accelerated expansion in a ACDM cosmology. The two conventions 
are simply related by f^^^ - [g(z = °°)/^(0)]/™^ where g(z) is the suppression factor defined as g(z) = D(z)(l + z) 
with D{z) being the linear growth function of density perturbations. In our numerical analysis of the CMB bispectrum 



in Section III as well as for the LSS analysis of Section IV we assume a flat ACDM cosmology and adopt the WMAP 
5-yr cosmological pai-ameter values: = 0.044,Q,. = 0.214,//o = 71.9,n, = 0.96,n, = 0,t = 0.087. For this 
choice one finds, for instance, /^f = 1.33/™'^. Since among our results are direct comparisons of CMB and LSS 
constraints, we will consistently assume the CMB convention for /nl throughout this paper 



B. Scale-dependence of the primordial bispectrum 

As implicitly assumed above, both for the local and equilateral type of non-Gaussianity, current observations only 
constrain the magnitude of /nl- However if /nl is large enough, it may be also possible in the near future to constrain 
its possible dependence on scale. There are well-motivated models of inflation, such as single-field models with a 
variable speed of sound, which naturally predict non-Gaussianity of the equilateral type with a scale-dependence of 
/nl 1361 l58l |62l. More recently, a model with local non-Gaussianity has been shown to also present a significant 
running of the amplitude of the primordial bispectrum, |48|. It can be well expected that for this kind of models, a 
measurement of, or simply a constraint on, this scale-dependence would provide additional limits on the fundamental 
parameters of the underlying high-energy theories. 

We remark, moreover, that allowing for a running of the non-Gaussian amplitude does not imply a degradation of 
the error on the amplitude itself. As we will discuss at length later on, for a given observable and a given form of the 
initial bispectrum, it is always possible to choose a proper pivot point to define the running parameter in such a way 
as to minimize, or indeed remove, any degeneracy between the two parameters. 

We consider a simple ansatz describing a mild scale-dependence of /nl- Specifically, we parametrize the initial 
curvature bispectrum for both local and equilateral models as 

/ \ "NG 

Bo(^i,^2,^3) = /nl^ FikukM, (II.6) 
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1^ \ "NG 

/nl^/nlI^ , (II.7) 



corresponding to the replacement 

f.„ fx,, I 
where we define 

K={kik2h)'^\ (II.8) 

and kp is a pivot point. The primordial bispectrum is therefore determined in terms of two parameters: the amplitude 
/nl and the new parameter «ng quantifying its running. 

We should point out that our definition departs from previous similar definitions introduced in the literature in two 
aspects. In the first, the overall scale K corresponding to the triangle formed by the wavenumbers k\, k2 and kj,, is given 
here by the geometric mean of the triplet characterizing the bispectrum configurations. Other works, e.g. lfT3ll27l . as- 
sumed the arithmetic mean, K - {k\ +k2 + ki)/3, which coincide with our definition only for equilateral configurations. 
This second definition in terms of a sum is indeed a more accurate and correct description of the bispectrum predicted 
by DBI models of inflation. Our choice is dictated by the much simpler numeri cal im plementation of the CMB es- 



timator that takes advantage of the separability of the geometrical mean in Eq. ( |II.8 1'. While we clearly recognize 
the importance of a full implementation of the proper form of bispectrum, which will be addressed in future work, 
we argue that we do not expect significant differences in our main results on equilateral non-Gaussianity due to this 
choice. We notice in the first place that the two definitions differ mostly for squeezed triangular configurations, while 
most of the signal is concentrated in equilateral ones, for the equilateral type of non-Gaussianity. We can quantify 
such diff'erence, following Babich et al. [49|, by introducing a generalized scalar product for the bispectra B, and Bj 
given by 

_ _ _^ Bi{kuk2,k,)Bj{kuk2,k,) 

lSP-B(ki,k2,h) 

where the sum runs over all triangular configurations given a specific volume and where represents the Gaussian 
variance of the bispectrum given by 

A'B(kuk2,k,) = -^P{ki)P(k,)P(k,), (11.10) 

with NT(ki,k2, k-i) being the number of fundamental triangular configurations in a given volume corresponding to the 
triplet ki, k2 and k^,. Notice that the variance A^B does not depend on the form of the bispectrum. In terms of the scalar 
product one can define the "cosine" 

B' ■ B ■ 

cos(/?„/?,)s , '' ^ , (11.11) 

V(B, ■ BdiBj ■ Bj) 

which can be interpreted as a measure of the correlation between the two forms. We assume a minimum value for k 
given by the fundamental frequency of a box the size of the Hubble volume, and a maximum value k^^^ - 1 /iMpc 
We compute the cosine between two bispectra both of the equilateral form, Fgq. and with a very large scale-dependence 
corresponding to hng = 1, but described respectively by the geometric and arithmetic mean in the definition of K. We 
find that the cosine between the two is larger than 0.99, indicating a strong correlation. In the case of the local form 
Floe, the correlation we find between the two definitions of K is the smaller value of about 0.91, a result expected 
since most of the signal comes from squeezed triangles where the difference between geometric and arithmetic mean 



For the Fisher matrix of large-scale structure observables, both forms can be equally considered and we will show that they both lead to similar 
outcomes in the equilateral case. 
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of the wavenumbers triplet is larger. In the local case, however, no explicit expression for a bispectrum with a strongly 
scale-dependent amplitude is given in the literature, so we assume our choice for to be a simple phenomenological 
ansatz, subject to improvements in future works. 

The second difference in our notation with respect to the literature regards the running parameter «ng- The notation 
hng was been first introduced by Chen [i36i and then adopted by Lo Verde et al. |13|, Chen et al. |58|, Khoury and 
Piazza l62l and also as nf^^^ in Byrnes et al. B8 l. to denote a running of the parameter /nl and it has been defined 
in analogy to the power spectrum spectral index as hng - 1 = dln\ fin^{k)\ /dink, where the scale-dependent /nl(^) 



could be defined, from Eq. (II.61 for equilateral configurations as /nl(^) = B(i,{k,k,k)/F(k,k,k). In the case of the 
spectral index, scale-invariance provides unity as the expected value for n^, in the case of a running non-Gaussian 
parameter there is no theoretical argument supporting an expected value for hng of one. In our definition a constant 
/nl corresponds to «ng = so that we can define, for equilateral configurations, 

«NG = — ^Tl • (11.12) 

dink 

We suggest this as a more natural definition to be adopted by future works on the subject. Lo Verde et al. ||T3l make 
use of the parameter k, the dimensionless time-variation of the speed of sound of inflaton perturbations^, related to our 
«NG by the expression hng - -2a-. 

Finally, we choose as the value for the pivot point kp = 0.04 Mpc ' , already adopted by ifTSll . Such choice is arbitrary 
but determines, given a specific probe and for a given model, the degree of degeneracy and correlation between the two 
non-Gaussian parameters. This choice indeed minimizes the degeneracy for the parameters of the equilateral model 
as determined by the bispectrum measurement in an ideal CMB experiment. We discuss the dependence of our results 
on the choice of kp in detail in Appendix [B| 



III. PRIMORDIAL NON-GAUSSIANITY IN THE COSMIC MICROWAVE BACKGROUND 

In this Section we describe the relation between the observed bispectrum of anisotropics in the CMB and the 
primordial curvature bispectrum, and we define the Fisher matrix for the two non-Gaussian parameters /nl and hng- 
We then show the expected results corresponding to the WMAP, Planck and CMBPol experiments together with those 
from the ideal CMB experiment. 



A. The Cosmic Microwave Background Bispectrum 

The harmonic coefficients of the CMB anisotropy a/„, = T ' J d^hAT(n)Y^^ can be related to the primordial fluctu- 
ation O as^ 



^<D(k)g^(^)7;„(k) + «^„, (III.13) 



where 0(k) is the primordial curvature perturbations, for a comoving wavevector k, gf(r) is the radiation transfer 
function where the index x refers to either temperature (T) or ii-polarization (E) of the CMB. A beam function bf and 



the harmonic coefficient of noise n^,,, are instrumental effects. Eq. (III. 13 1 is written for a flat background, but can 
easily be generalized. 



^ We denote this quantity witli tlie symbol s in Section VI 

^ Note tliat the formulae in this Section are written in the opposite Fourier Transform (FT) convention than the one used for formula jll. 1 1 and that 
will be used in all the other Sections of this paper. The change of convention for this Section was adopted in order to be consistent with previous 
literature in the non-Gaussian CMB field. 
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Any non-Gaussianity present in the primordial perturbations C>(k) gets transferred to the observed CMB via 



Eq. (III. 13 I. The most common way to look for non-Gaussianity in the CMB is to study the three-point function 
of temperature and polarization anisotropics in harmonic space. Such quantity is called the CMB angular bispectrum 
and is defined as 

fif;, ={a'l 4 a'), (III. 14) 

{lt213,'"l'K2'"3 'inil Cl«il tint}' ' 

and the angular-averaged bispectrum is 



nPir 



/—i \ OTi TO2 '"a ' (i(2h,m,m2m3 ' ^ ' 



which using Eq.( |III. 13] l can be written as: 

rpi' - (A-n-\^( ;\e>+e^+e^ V / ^3 \ f dh^d%^cPk^ ^ ^ ^ 

BtAi. - L mi m2 m,] \ (2nf (2nf (2;r)3 ^ tunS^^^^ hmS^^)^ e,,nS^i) 

Xgi (ki)gl{k2)gl{h) <(D(ki)0(k2)0(k3)> , (III.16) 



where (<l)(ki)<I)(k2)<l>(k3)) is the primordial curvature three-point function as defined in Eq. (II. 1 1. 

To forecast constraints on scale-dependent non-Gaussianity using CMB data, we will perform a Fisher matrix anal- 
ysis. We will consider some non-zero /nl as our fiducial value for the Fisher matrix evaluation. Clearly, in order to 
be able to constrain a scale-dependence of /nl, its amplitude must be large enough to produce a detection. If /nl is 
too small to be detected (/nl < 2 is a lowest theoretical limit even for the ideal experiment), we will obviously not be 
able to measure any of its features, either. In the following we will then always consider a fiducial value of /nl large 
enough to enable a detection. Following IT] |2l [63J, the Fisher matrix for the parameters /nl and hng (for generic 
fiducial values /nl 5^ and «ng ^ 0) can be written as: 

Z Z ^^(cov").. (ni.i7) 

The indices a and b can take the values 1 and 2 corresponding to the Fisher matrix parameters p\ = /nl and p2 = «ng, 
respectively. Indices ijk and pqr run over all the eight possible ordered combinations of temperature and polarization 
given by TTT, TTE, TET, ETT, TEE, ETE, EET and EEE; the combinatorial factor A(^(,j^ equals 1 when all ^'s 
are different, 6 when {\ - li- (3, and 2 otherwise. The covariance matrix Gov is obtained in terms of Cj^, Cf^, and 
Cj^ (see imiU) by applying Wick's theorem. 



For non-Gaussianity of the equilateral type, for which the functional form F(ki,k2,ki) is given by Eq. (II. 5 1, we 
have 

, ijk 



^^tkh ^ / (2A + 1)(2^2 + mh + l) l£i h h 
5/nl V 4;r ^0 

X 6 J" 'P-dr [-a[fi{fi\ + 2 perm. + j6'f,ri<Jl + 5 perm. - Id^d^/f^] , (III.18) 



where the functions a,j6, 6, and y are given by: 
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2 r Ik \"'^°^^ 

- ykeg',{k) jt(kr)y-^ , (III. 19) 

2 r Ik Y^'^l^ 

-ykk~'g\{k)jt{kr)y-^ AoT'-i, (III.20) 

2 r / yfc \"'^"^^-' 

y[{r) = - ]^dkkg[{k)ji(kr)y-'^ A'^H'^"'-'''', (III.21) 

2 r Ik \"'^°^^ 

6',(r) ^ -jdkg^(k)Mkr)\^-j Af fc2(«,-i)/3_ (^ 22) 

In the expression above we use the dimensionless power spectrum amplitude Ao, which is defined by f o(^) = 
A(p^"^^*"'"'\ where n., is the tilt of the primordial power spectrum. 

The derivative with respect to «ng produces an extra In(fe) factor in the integrals over the CMB transfer function. 
The explicit form is given by 

^^e!L /(2^, + 1)(2^2 + 1)(2^3 + 1) / ^1 {2 {3 



4(r) 
Air) 



driNG V 4-n \ 

/-*CC 

*J 

AyiA +P'cMA +f^tjiA + ^ Pem. - 2~6',/,/,^ - 25[6{/,^ - 16\6{5\\ , (III.23) 
comprising a total of 30 additive contributions, where the functions a,^, S, and -y are given by: 

2 r^°° I k Y"^^^^ I k \ 

a-f(r) = — j^^ dkk^g;(k)je(kr)\^-j InJ-j (III.24) 

\nNG/3 



2 r^°° I k \ I k\ 

p-^(r) = — j^^ dk k-' g',(k) je(kr)\^-j InJ-jAoT-' (III.25) 

2 r^°° / k \"'~'°'-' I k\ 

-I dkkg',ik)Mkr)[-] ln(^)A^/^^.<-»/^ (III.26) 



m =- |^p..[(.);X..)(^f"^ln(A)A-.---. (III.27) 

In a similar way, from Eq. ( II.4 1, one can derive the following expressions for the bispectrum derivatives in the local 
case. 



^^VA / (2A + 1)(2^2 + 1)(2^3 + 1) / A {2 r ioJok 1 ,TTT'78^ 

= y ^ ^ Oj^l'-dr [-o^eAPe, + 2 pem.J , (III.28) 



and 



^^'eih ^ I (2^1 + 1)(2^2 + 1)(2^3 + 1) / ^1 {2 £3 
5«NG V 47T ^ 

r^flfr [4,y3l/??3 + a'f^/,^ + a^fii}',^ + 2 perm.] . (III.29) 

We will not consider a marginalization over cosmological parameters. For a discussion of the effect of uncertainties 
on the cosmological parameters on the determination of the non-Gaussian amplitude parameter /nl from CMB mea- 
surements, see 1*641. We modified the publicly available CMB fast code 1651 to compute the transfer functions and 
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B. Constraints from the CMB Bispectrum 

We compute here the expected uncertainties on the two non-Gaussian parameters /nl and «ng from the Fisher 
matrix analysis of the CMB bispectrum assuming as fiducial values /nl = 50 and hng = for the local model and 
/nl - 100 and hng = for the equilateral one. We discuss in detail the dependence of these results on the choice of the 
fiducial values in Section [riI C | We consider the specifications for the current Wilkinson Microwave Anisotropy Probe 
[WMAP. I661 and Planck |67 1 missions, the proposed satellite mission CMBPol |68 69 1 and an ideal CMB experiment, 
assuming a full-sky coverage. We consider the CMB bispectrum up to /"^ax = 1200 for WMAP, ^max = 2500 for Planck 
and ^max = 3000 for CMBPol and the ideal case. 

In Fig. [T| we compare the expected 1-cr contours on local model parameters /^°j^ and hng for WMAP (upper left 
panel), Planck [upper right), CMBPol (bottom left) and for an ideal CMB experiment (bottom right). Dashed blue 
lines coiTespond to the limits from the temperature information alone, dotted green lines to polarization (BEE), while 
the continuous red lines correspond to all bispectrum combinations. The figure assumes for the pivot point kp = 0.04 
Mpc resulting in a significant degeneracy between the two parameters. As shown in Appendix [b] a better value for 
the local model CMB bispectrum would have been closer to kp - 0.02 Mpc ' for Planck and kp = 0.03 Mpc ' for 
CMBPol, indicating that most of the signal is coming from scales larger than the pivot point. An optimal approach 
would actually require one to determine the pivot scale not only depending on the shape of non-Gaussianities, but also 
according to the specific experiment and its multipole range and noise properties. However, it will be shown later that 
by combining CMB and large-scale structure observables the degeneracy between riNc and /nl is largely removed. 
Moreover, the pivot scale kp - 0.04 was the one adopted in [ 13], and it thus allows an easier comparison between the 
two results. For these reasons we will always assume kp - 0.04 Mpc"' throughout the rest of this work. In Fig.0 
we plot the expected marginalized (thick, blue lines) and unmarginalized (thin, red lines) 1-cr errors on the individual 
parameters /nl (left panel) and mng (right panel) for Planck (upper panel) and CMBPol (lower panels) as a function 

of l/nax- 

We find that WMAP and Planck should be able to provide a 1-cr uncertainty on the running of /Jjj^ respectively 
equal to 

A«NG ^ 0.68 4^ — [WMAP], (III.30) 

/nl yfsky 

and 

A«NG ^ 0.10 4^ — [Planck], (III.31) 

/nl V/*' 

after marginalizing over /nl, while for CMBPol we find a smaller error by about a factor of two 

Ahng ^ 0.05 4^ — ^ [CMBPol]. (III.32) 

/nl V/*' 

Notice that the marginalized error on hng does not depend strongly on the choice for the pivot point kp. The un- 
marginalized 1-cr errors, for a fiducial f^- = 50 and f,ky = 1, are, in fact, Ahng — 0.24, 0.07 and 0.04 for WMAP, 
Planck and CMBPol, respectively. 

To give a sense of the magnitude of the errors provided by the Fisher matrix analysis, we can consider the change in 
the non-Gaussianity amplitude /nl(^) over the range of scales probed by the CMB (which roughly goes from 0.001 to 
0.1 /;"' Mpc), and compare it to the /nl sensitivity of a given experiment. We roughly expect a departure from nfjc = 
to be undetectable if the variation it produces in /nl(^) over the available range of k is below the /nl uncertainty for 
the experiment under examination. For example, taking a fiducial value /nl = 50 at = kp, a running of hng - 0.1, 
corresponds to the two extremal values /nl(^ = 0.001 /iMpc"') ^ 33 and /nl(^ = 0.1 /iMpc"') ^ 53. At smaller 
scales, such as those probed by the cluster abundance, we have /nl(^ - 0.5/iMpc ') ^ 62. These variations are well 
within current uncertainties on the amplitude parameter Our previous argument then suggests that hng = 0.1 is too 
small to be detected by WMAP, in agreement with the full Fisher matrix analysis presented above. For a larger value 



10 




FIG. 1: Local model. 1-cr constraints on /nl and wng assuming kp = 0.04 Mpc"' and fiducial values /nl = 50, ;jng = 0. Dashed 
lines correspond to the limits from the temperature information alone, dotted lines to polarization (EEE), while the continuous lines 
correspond to all bispectrum combinations. We consider WMAP (upper left panel), Planck {upper right), CMBPol {bottom left) 
and an ideal CMB experiment {bottom right). 



of the running parameter as hng 
while for the negative value hng - 



0.4, one obtains ft^i^(k = 0.001 /; Mpc"') 10 and /nl(^ 
-0.5, /NL('t = 0.001 h Mpc"') ^ 250 and ft^idk 



0.1/iMpc"') 



O.l/iMpc"') ^ 63 
^ 40. This amount 



of running starts to be in the detectability range of WMAP, again in agreement with formula of Eq. ( 111.30 1. We notice 
that the large value obtained for the largest scale under the assumption of a running close to the expected WMAP error 
on hng, is close to the uncertainty derived in |5 1 for the smallest C-hin in their analysis of the WMAP bispectrum. 
Fig. |3]shows the same results as Fig.[T]but for the equilateral model. The reduced degeneracy with respect to the 
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FIG. 2: Local model. Expected marginalized (thick, blue lines) and unmarginalized (thin, red lines) errors for /nl (left panels) and 
wng (right panels) as a function of l„,ax for Planck (central panel) and CMBPol (lower panels). 



local case indicates that the chosen pivot kp = 0.04 Mpc"' is close to the optimal value for the equilateral model. This 
also reflects the fact that the overall signal is, for this model, relatively larger at smaller scales. This is due to the 
statistical weight of the number of configurations close to equilateral which increases faster with l,„ax then the number 
of squeezed configurations correlating small scales to the few available large scales, relevant for the local bispectrum. 
In Fig.|4]we plot the expected marginalized (thick lines) and unmarginalized (thin lines) l-cr errors on the individual 
parameters /nl (left panel) and hng (right panel) for Planck (upper panels) and CMBPol (bottom panels) as a function 
of Imax in the equilateral case. 

The 1 -cr uncertainties on the running parameter «ng expected from WMAP, Planck and CMBPol are given by 

100 1 

AnNG ^ 1-1 — —j= [WMAP], (III.33) 
AnNG ^ 0.30 ^ — ^ [Planck], (III.34) 




A«NG ^ 0.17 — [CMBPol]. (III.35) 

/n 

The corresponding unmarginahzed values assuming - 100 and fsky = 1 are A«ng - 0.47, 0.28 and 0.16 for 
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FIG. 3: Equilateral model. 1-cr constraints on /nl and hng assuming k,, = 0.04 Mpc"' and fiducial values /nl = 100, hng = 
0. Dashed lines correspond to the limits from the temperature information alone, dotted lines to polarization (EEE), while the 
continuous lines correspond to all bispectrum combinations. We consider WMAP (upper left panel), Planck {upper right), CMBPol 
{bottom left) and an ideal CMB experiment {bottom right). 



WMAP, Planck and CMBPol respectively. In the equilateral case, the larger uncertainties on hng allow for a larger 
variation of the amplitude of non-Gaussianity over the probed range of scale. However, we would like to remind the 
reader that the larger uncertainties in /nl, "ng for the equilateral case are only due the specific normalization choice 
for the equilateral bispectrum, and do not reflect a larger detection power for local non-Gaussianity with respect to 
equilateral non-Gaussianity |50|. Assuming «ng - 0.3 and fm^{k - hp) - 100, we have /nl(^ = 0.001 /iMpc"') ^ 30 
and ftii^ik = 0.1 /z Mpc"') ^ 120 while for «ng = -0.3, /nlC^ = 0.001 /iMpc"') ^ 330 and f^^k = 0.1 /i Mpc"') ^ 84. 
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'max 'max 

FIG. 4: Equilateral model. Expected marginalized (thick, blue lines) and unmarginalized (thin, red lines) errors for /nl C^/j panels) 
and wng (right panels) as a function of l,nax for Planck (central panel) and CMBPol (lower panels). 



These are, again, relatively small variations, stressing the constraining power of the CMB bispectrum on the running 
parameter. 



C. Dependence on the fiducial values 

The dependence of the previous results, both for the local and equilateral models, on the fiducial value f^^ of /nl 
is trivial since for the Fisher matrix, Fij, with /, j - I corresponding to /nl and 7 = 2 to n^c, we have the scaling 
F\2 ~ /nl and F22 ~ /^l- While the errors on /sjl clearly do not depend, in our approximation, on the assumed value 
for the same parameter, for the errors, both marginalized and not, on n^c we have in general 

A«NG ~ 1//nl' an-36) 



as emphasized by the expressions in Eq.s (III.30 1 to (III.35 1. The correlation coefficient between the two parameters, 
defined as C12 = {F^^)i2/ ')ii(-P'"')22 is independent of f^^^ so that the value of the latter does not affect the level 
of degeneracy. On the other hand, the fiducial value n^^ of hng clearly has a large impact on the distribution of the 
signal over the relevant range of scales. We can therefore expect it to determine in part the level and direction of the 
degeneracy between /nl and «ng- 
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FIG. 5: l-cr contours for CMBPol assuming different fiducial values for the parameters. Leftpanel: local model for /nl = 5,25,50. 
Right panel: equilateral model for /nl = 10, 50, 100. In both cases we also consider n^Q = (continuous, red lines), 0.1 (dashed, 
blue lines) and -0.1 (dotted, green lines). 



These considerations are illustrated in Fig. |5] where we compare the contours for temperature and polarization 
shown in Fig.[T]and Fig. |3]to the same contours derived assuming dilferent fiducial values for /nl and «ng. On the 
left panel we show the local model results for /nl = 5, 25 and 50 and for hng = and +0.1. On the right panel we 
show instead the equilateral model for /nl - 10, 50 and 100 with «ng - and +0.1. In this last case, the effect on the 
parameter degeneracy is particularly evident, although we would like to stress again that, once a given experimental 
configuration has been chosen, the degeneracy can always be removed through an appropriate choice of pivot and it is 
not going to be an issue in the experimental determination of /nl and hng- 



IV. PRIMORDIAL NON-GAUSSIANITY IN THE LARGE SCALE STRUCTURE 

The impact of non-Gaussian initial conditions on the evolution of the matter and galaxy distributions and their effect 
on large-scale structure observables has been the subject of several works over the last two decades focusing, for the 
most part, on the galaxy clusters abundance |6, 7, 8 9 10, 70 1, on higher-order moments of the density probability 
distribution lfT5l[T6HT7ll . on higher-order galaxy correlation functions, IITSl [191 12011211 and the Minkowski Functionals 
of the density field ll30l |3TI . More recently, a strong effect due to non-Gaussian initial conditions of the local type 
on the bias of dark matter halos has been measured in A^-body simulations by Dalai et al. [22J. This result, somehow 
anticipated by earlier works on high-peak clustering |71 72 1, spurred a series of works Il23l l24l [25l l26l [27l [28l 1611 
[73] |74||75]|76| indicating that significant constraints on the amplitude of local non-Gaussianity, comparable to those 
obtained from the CMB bispectrum, can be derived from measurements of the power spectrum of galaxies and quasars 
in current large-scale surveys |2T 26]. At the same time, further studies of the effect of non-Gaussianities on the halo 
mass function are being proposed |77, 78„ 79. 80|. 

We will not attempt here to perform a complete analysis of LSS probes including all relevant observables (and their 
covariance) as this would be beyond the scope of this work. We will derive, instead, simple estimates of the expected 
errors on the amplitude and running parameters /sjl and hng, choosing specifically a representative observable for 
each model: namely the bias of galaxies as determined from power spectrum measurements in galaxy surveys for the 
local non-Gaussian model, and the galaxy bispectrum for the equilateral model. As we will argue in the next Section, 
we expect the two different models to affect in quite a different way the large-scale clustering of biased objects. We 
therefore make, in choosing these specific and different measurements, what we consider a reasonable assumption on 
the impact of different non-Gaussian models, warning that further theoretical work, possibly supported by numerical 
simulations, will be needed to strengthen our understanding of this phenomenon. 
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Regarding the possibility of a scale-dependent non-Gaussianity, Lo Verde et al. ifTSll provided the first analysis 
of the possibility of constraining a running /nl parameter by combining current limits from the CMB with future 
measurements of cluster abundance. Focusing in particular on the equilateral model for the curvature bispectrum, 
this work assumes the amplitude of /nl(^) = /nlC^/^p)^* to be constrained by the CMB bispectrum at the pivot 
point kp - 0.04 Mpc ' and derives the expected constraints on the running parameter"* s by considering the effective 
amplitude of /nl(^) at the smaller scales (k ~ 0.3 - 0.6 /z Mpc"') probed by cluster surveys (see Fig. 1 in fU l). For an 
all-sky cluster survey up to redshift Zmax = 1-3 with single mass bin defined by the threshold M > 1.45 x lO'"* /i"' Mq 
they find, in terms of our running parameter «ng = -2i, the l-cr constraints, marginalized over Q„,, erg and h, assuming 
the fiducial values /nl - 38 and /ing - 0, Ahng - 3.38 with a WMAP prior A/nl(^ = kp) - 150, and A«ng - 2 with a 
Planck prior Afyi^(k - kp) = 40. Note that their analysis does not include the simultaneous limits that measurement of 
the CMB bispectrum alone is expected to provide on both the amplitude /nl and running hng that have been discussed 
in the previous Section. In terms of these fiducial values, Eqs. ( ni.33| l and ( 111.34 1 give, for instance, Amng - 2.9 



for WMAP and Ahng - 0.8 for Planck, a factor of 1.2 and 2, respectively, better than the results of fT3l. It is to be 
remarked, however, that in the event of a detection, an independent confirmation from an observable probing quite 
diff'erent scales such as the cluster abundance would be of great importance. 

The matter overdensity in Fourier space is related to the curvature perturbations Ok at early times during matter 
domination by the Poisson equation as 

(5k(z) = M(^,z)(5k, (1V.37) 

where we introduced the function 

with T(k) being the matter transfer function and D{z) the growth function. The linear matter power spectrum is 
therefore given by 

Poik)^ M^{k,z)P^{k). (IV.39) 

A departure from Gaussianity in the initial conditions results in non-vanishing higher-order, connected, correlation 
functions for the matter field proportional to the corresponding correlation functions for the curvature fluctuations. In 
particular, for the initial matter bispectrum Bq we have 

BQ{kuk2,h) = M{kuz)M{k2,z)Miki,z)B<i,{kuk2,h). (IV.40) 

For simplicity, we will ignore the non-Gaussianity represented by a non-vanishing initial trispectrum — and by higher- 
order correlation functions — although its effect on the galaxy bispectrum can be significant, see 1371 138 1. On the other 
hand, their expressions involve extra parameters, other than /nl, and their scale-dependence deserves a specific study. 

In the rest of this Section we will summarize and discuss in some details the recent literature regarding the effect 
of non-Gaussian initial conditions on the bias of galaxies and clusters and we motive our choice for the observables 
considered for our simple analysis. We will then derive the expected uncertainties on the two parameters for the two 
models of the curvature bispectrum introduced in Section[lI]considering measurements of the galaxy power spectrum 
for the local model and the galaxy bispectrum for the equilateral model. 



A. Halo bias and primordial non-Gaussianity 

As first shown by Dalai et al. Il22l and later confirmed by Grossi et al. |29|, Desjacques et al. [74|, Pillepich et al. 
||75l in numerical simulations with non-Gaussian initial conditions of the local kind, the large-scale bias of halos can be 



The dimensionless variation of the speed of sound s was denoted k in 1131 . 
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greatly affected by relatively small values of /nl- Dalai et al. ^T\, Slosar et al. Il24l . Afshordi and Tolley ll26ll provided 
expressions for the halo bias in the presence of local non-Gaussianity based on the peak-background split formalism 
IISTl . They find the correction Abh, i (k, M, z, /nl) to the Gaussian linear halo bias bh, i (M), defined as b'^'^ = bh, i + Abh, i , 
to be given by 

Abh,iik,M,z,fNL) = 2fnL[bh,i(M,z) - m ^Jr-, (IV.41) 

M{k, z) 

where 6c - 1.68 is the critical overdensity for spherical collapse and /nl - 1.33/nl is the non-Gaussian parameter in 
the large-scale structure convention. The strong scale- and redshift-dependence of such correction allowed Slosar et al. 
Il24l and Afshordi and Tolley 1261 to place significant constraints on f^^- already from present large-scale structure 
data. 

A different approach has been followed instead by Matarrese and Verde 1231 and by Taruya et al. |27| to derive 
a similar expression for the correction to the power spectrum of biased tracers of the dark matter distribution. In 
particular, Matarrese and Verde f23l considered a peak number density defined by the expression Py,,M(x) = 0[6r(x) - 
6c], where is a step function. Based on earlier works by Grinstein and Wise |71| and Matarrese et al. |72|, they 
derive the correction to the peak correlation function in terms of the three-point correlation function. Their result can 
be rewritten in the form of a linear peak bias correction given by 

Ab„^i(k, M, z, /nl) = 4 ^TiAt^ r d'qBodk, q, |k - q|) (1V42) 

where v = 6c/<Tr represents the peak height with ctk being the r.m.s. of matter fluctuations filtered on the scale R 
and where Po,r - M^ik, z)P<s,(k) and Bo,R(fci, ^2, ^3) - MR{ki,z)MR{k2,z)MR{ki,z)B<s,{ki,k2,k{) are the initial, filtered 
matter power spectrum and bispectrum and where we introduce, as short-hand notation, MR{k, z) = WR{k)M{k, z) with 
WR{k) being the smoothing function. In the large-scale limit (k — > 0) and assuming the local form for the curvature 
bispectrum Bij, with no running, one can use the approximation 



r d'qBo,Rik, 9, |k - q|) - 4/nl , (IV43) 

\,R(k) J MR{k,z) 



\ d^aRn Jk. Ik - nl^ 4 — - 

Pom 

which leads to the correction to the linear peak bias 



Abp,dk,M,z,fnL) - 2/NL477Tr^ = 2/nl(^,i - D^.tttt^' aV.44) 
criMR(k,z) MR(k,z) 



equivalent to Eq. (|1V.41 1 after identifying bp^i = 1 -H v/cr« as the Eulerian linear peak bias. 



Taruya et al. |27 1 computed the 1-loop corrections to the galaxy power spectrum due to a non-vanishing primordial 
bispectrum, assuming a local bias relation between galaxies and the smoothed matter density 6r(x), 

6gix) = bi6R(x) + ^b26l(x) + ... , (IV45) 
finding a correction to the linear bias parameter that can be written as 

Abi(k, M, z, /nl) = b2 r d'qBoMk, q, |k - q|) aV.46) 
which, in the large-scale limit becomes 

Abi(k, M, z, /nl) - 2/nl^2 -tttt—, ' (1^.47) 
MR(k, z) 



and which notably involves the quadratic bias parameter b2 from Eq. (1V.45 1. 
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The three expressions of Eq. ( |IV.41| l, Eq. ( |IV.42[ ) and Eq. ( |IV.46| l all show the same scale- and redshift-dependence 
but with a different constant of proportionality which depends on the specific assumptions. In other words, references 
Il23l and (|27 1 derive an effect on the halo/galaxy bias that, for the local model, can be shown to be equivalent, at large 
scales, to the one described by Eq. (1V.41 1, and they actually coincide if one considers as well the high-peak limit 
v/cr^ » 1 since, in this case 



[bh,i - l]6c - b2crl ^ = -|, 



(IV.48) 



once one assumes an expression such as the one derived by Mo and White 11821 for the galaxy quadratic bias parameter. 
On the other hand, both |23| and ETll provide an expression that can be applied to any model of primordial non- 
Gaussianity as they depend generically on the primordial curvature bispectrum. For instance, as pointed out in |27i. 



while for a local non-Gaussianity both Eq. ( 1V.42 1 and Eq. ( 1V.46 1 predict a large effect as the one measured in A^-body 
simulations, in the equilateral case the effect is largely suppressed. As shown in Fig. 6 of ll27l . the relative correction 
to the Gaussian bias parameter at redshift z = 1, with b\ - 2 and b2 - 0.5 stays below 1%, even when a running of 
/nl corresponding to «ng - 0.6 is considered. 

It should be noted that while the parameterization of a scale-dependent /nl(^i , ^2, ^3) described in Sectionjll] which 
is to be understood as a parameterization of the initial bispectrum, can be easily and unambiguously implemented in 
Eq. (1V.42| and Eq. (lV.46i, the same is not true for Eq. (IV.41 1. In fact while Eq. (1V.41 1 suggests one replace a 
constant /ml as indicated by Eq. ( 11.7 1 with k being the unequivocal scale of the problem, in the other case, assuming 
k - (A;i^2^3)'^"^ in the expression for the curvature bispectrum, the scale dependence is partially integrated over, 
leading, as we will show in detail in the next Section, to a quite different behavior. Based on these considerations we 
assume for our analysis an expression for the correction to the linear bias parameter proportional to the integral of 



Eq. (IV. 46 1 and Eq. (IV. 42 1: this provides an expression for the correction to the galaxy power spectrum in the case of 



local non-Gaussianity and, at the same time, a plausible argument for neglecting altogether any effect on the galaxy 
power spectrum in the equilateral case. 

For equilateral non-Gaussianity we consider instead its effect on the galaxy bispectrum, and, in particular, we assume 
such effect to be given simply by the primordial component, neglecting corrections due to non-linear galaxy bias. This 
is a quite conservative choice. In fact, it has been recently shown by Jeong and Komatsu |37| and by Sefusatti ll38ll 
that significant corrections at large scales, similar to those affecting the galaxy power spectrum, are expected as well 
for the galaxy bispectrum. Since these early results still require a comparison with numerical simulation to properly 
test their theoretical predictions, we will consider for equilateral non-Gaussianity the sole initial component present as 
well in the matter bispectrum. 

A complete analysis of large-scale structure data should naturally include two-point statistics as well as higher- 
order correlators. In this work, however, we neglect for the time being the extra information that could be extracted, 
for instance, from the galaxy bispectrum on local non-Gaussianity, leaving this subject for a future work. As for 
the equilateral model, we consider as a large-scale structure probe the galaxy bispectrum, extending the analysis of 
Sefusatti and Komatsu lISTI to the case of a running f^, keeping in mind that this analysis needs to be improved 
by taking into account the large-scale corrections due to non-local bias studied in If37l l38l that will probably lead to 
smaller uncertainties on the non-Gaussian parameter /nl both for the local and equilateral models. 



B. Constraints on local non-Gaussianity from the galaxy power spectrum 

In this Section we perform a Fisher matrix analysis to obtain an estimate of the constraints achievable from power 
spectrum measurements in future galaxy surveys on and «ng- We assume the correction to the galaxy linear bias 
parameter b\ to be given by 

A/7i(^,z,/nl) ^2q f^dk) (bi - 1) S,-^— ] (IV.49) 

Mnik, z) I + (ak)" 

This expression assumes the normalization of the bias correction to be given by comparison to A^-body simulation and 
encodes the actual dependence on the scale-dependent curvature bispectrum in the function fyhik)- Specifically, we 
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include, to the theoretical expression of Eq. ( |1V.41| |, the correcting factor q - 0.75 introduced, and motivated in terms 
of the ellipsoidal collapse model, by Grossi et al. \29\ and the small-scale suppression factor + {ak)"] derived by 
Desjacques et al. 11741 from a fit to numerical simulations at smaller scales, with a - 23.5 and n - 1.69 which we will 
assume valid for all redshift and densities. Such damping has been observed as well by Pillepich et al. ifTSl . which 
also provided a fitting formula consistent with the one we assume on relevant scales. In other terms, we are assuming 
the amplitude of the bias correction to be proportional to [bi{M,z) - 1]/[1 + {ak)"} since this factor has been compared 
to simulations in the case of a constant /nl. However, we assume at the same time a scale-dependence consistent with 
Eq. (|1V.46|), so that 



/nl(^) = 



_MR(k,z) 
4o-'- 



iPoAk) J 



d^qBo.R(k,q,\^-q\) 



4^1 I 



d^qMR(q)MR{\k - q|) 



B^{k, q, |k - q|) 



d'qMR(q)MR(\k - q|) 
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NL 



'(kq\k-q\) 



1/3 



P<!,(q) + Po(|k - q|) + 



Po(g)f o(|k - q|) 
P<s>(k) 



av.50) 



where the last equality assumes the curvature bispectrum of the local model. The dependence on the smoothing scale 
R is very weak at large scales if «ng - (see also Fig. 1 in fSj]), since in this case, for the local model, we have 



/nl 
2crl 
/nl 
2cr' 



j d'qMRiq)MR(\k-q\) 

\ [ d\Ml{q)2P<s,{q) = /i 
R J 



P^(q) + P<!,(\k-q\) + 



Po(?)Po(|k - q|) 



P<bik) 



(IV.51) 



but it is significant, even at large scales, when the running parameter is different from zero. In Fig. |6] we plot the 



scale-dependent part of the bias correction, fy]^{k)/M(k,z) according to Eq. (|1V.50|) assuming the amplitude f^^- - 1, 



the pivot kp - 0.04 Mpc"' and a smoothing scale R = 5 Mpc (left panel) and R = 10 Mpc (right panel). Thick 
lines show the behaviour for the different values «ng = (continuous line), hng = +0.5 (dashed) and hng = -0.5 
(dotted). The thin, dashed vertical line indicates the position of the pivot point, kp in h Mpc 

In the left panel of Fig. |6] the thin (gray) lines show the behaviour corresponding to a simple substitution of /nl 
with fNL(k) = fNL(k/kp)"'"^ in Eq. (1V.41 1. One can clearly see that such naive substitution would lead to a very 
different effect on halo bias, particularly on the relevant range of scales. The integration in Eq. ( |1V.50| | in general tends 
to dampen the effect of a running /nl since most of the scale-dependence is integrated over Moreover, it induces an 
"effective" pivot scale that depends on the value of R and it is greater, in momentum space, for larger values of R. 
This can be understood in terms of the expected relative value of the integral in Eq. ( 1V.50 1, at the pivot point kp with 
respect to the value of ctr. Since both the chosen values R - 5 /; ' Mpc and R - 10/i"' Mpc correspond to momenta 
much larger than kp = 0.04 Mpc"', for the smaller value of R - 5 h^^ Mpc we expect the integral to provide a larger 
value for positive hng relative to ctr than the case of R = 10 /i ' Mpc. The dependence on the assumed value of kp 
is instead more trivial, as a shift in kp corresponds to a shift in the effective pivot point of the bias correction. It is 
interesting to note that this peculiar behaviour might be relevant to correctly identify a theoretical description of this 
effect, once the scale-dependence of the curvature bispectrum is properly implemented in numerical simulations. 

In our analysis we will assume for simplicity a single representative value ofR = 5 /z"' Mpc, although the smoothing 
scale should in principle be related to the halo population of interest. We notice, however, that while the degeneracy 
between the two non-Gaussian parameters is affected to some extent by this choice, the overall value of their marginal- 
ized uncertainties does not change significantly for our choice of the pivot scale. 

For a given survey, defined uniquely by its redshift range, sky coverage and expected galaxy number density iig, 
we consider shells in redshift of maximum size Az - 0.5 and evaluate all quantities at the mean redshift zj- This 
implies that, even in the ideal case of an all-sky survey, the volume of each shell does not exceed 100 h^^ Gpc^ for our 
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FIG. 6: Scale-dependent factor, fyi^(k)/M(k,z), of the bias correction for /nl = 1 and hng = (continuous red line), «ng = 0.5 
(dashed red line) and «ng = - 0.5 (do tted red line). The vertical dashed line indicates the pivot scale assumed, kp = 0.04 Mpc"'. 
Left: smoothing scale of Eq. ilV.SOi, R = 5/i"' Mpc. Right: R = 10 h^' Mpc. The thin gray lines on the left panel show the 
behaviour correspondent to a simple substitution of /ml with /nl(^) = (k/kp)"'^'^ in Eq. l |IV.41^ . 



cosmology, so that the largest scale probed corresponds to a wavenumber k,„i„ ~ 0.0013 Mpc~^ We conservatively 
ignore large-scale correlations among dilFerent redshift bins. 

The Fisher matrix for the non-Gaussian parameters and is therefore given by 

dPJki,Zi)dPJki,Zi) 1 

where A^^ is the number of redshift bins and A^^ is the total number of wavevectors kj from kf - In/V^^^ to a maximum 
kmax - 0.03 h Mpc"' in steps of ^/ in the given redshift shell. 

We will consider, again conservatively, only linear corrections in /nl to the galaxy Gaussian linear bias bi according 



to Eq. ( |IV.49| l, so that 

Pg(k) ^b][l+ 2Abi(buk,z,hL)] P(k) av.53) 

with P{k) being the linear matter power spectrum. The power spectrum variance AP^{k) is approximated by its Gaus- 
sian component 



,2 



^^AAbi{kJ^i}P{k) 



PUk) 

where f ,o, is the Gaussian galaxy power spectrum including shot-noise 



(IV.54) 



1 



PUk) = Pgik) + (IV.55) 



with Hg representing the galaxy density of the survey. The last factor in Eq. (IV.54 1 corresponds to the linear correction 
to the Gaussian variance due to effect on bias of a non-vanishing /nl, which we will assume as fiducial value. We 
ignore for simplicity corrections due to higher-order correlations. The galaxy bias parameter is assumed to be known 
and it is computed by means of the halo model prescription and halo occupation distribution described in l,21J . to 
which we refer the reader, which in turn is based on ll83l l84l [85l [86l ISTl l . 

In Table|l]we show the specifications of LSS surveys that we will assume for the analysis in the following Sections. 
For the power spectrum analysis we consider both spectroscopic and photometric surveys, under the assumption 
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that photometric uncertainties have a negligible impact on the determination of galaxies clustering properties on the 
interesting range of scales. Among the surveys, proposed or under construction, we choose those able to provide 
constraints on the non-Gaussian parameters, comparable to those provided by future CMB observations. 

In Fig.|7]we show the 1-cr contours of the uncertainty on the non-Gaussian parameters and hng determined from 
the Fisher matrix analysis of the galaxy power spectrum assuming fiducial values = 50 and hng = 0. On the left 
panel we consider an ideal all-sky survey with a galaxy number density of 10"^ Mpc""' up to redshift 3. We plot the 
contours corresponding separately to the redshift intervals < z < 1 {dotted line), 1 < z < 2 (dot-dashed), 2 < z < 3 
(dashed) and < z < 3 (continuous). On the right panel, we show the contours for a few proposed large-scale surveys. 
The unmarginalized errors on the amplitude parameter /nl are larger by a factor of a few with respect to the results of 
the similar analysis of lISTl . depending on the survey characteristic and specifically on the value of the linear Gaussian 
bias bi. The reason for this is maily the inclusion in the variance for the power spectrum of the correction to the bias, 
since we are assuming a fiducial value for /nl different from zero. This reduces the signal-to-noise at the largest scale. 
We notice that, for a fiducial /nl - 0, we are able to recover the results of |61| only replacing the factor (bi - 1) in 
Eq. (IV.49 1 with bi, as it seems to be assumed by Eq. (9) defining the Fisher matrix in |61 J. 

The direction of the degeneracy displayed by all contours reflects the choice for the pivot point kp and indicates that 
most of the signal is coming from scales larger than the pivot scale. In fact, we remind the reader that the Fisher matrix 
analysis is limited to 0.001 Mpc < A: < 0.03 /i"' Mpc by choice and this results in the quite uniform results for the 
diff'erent surveys in the left panel of Fig. [7] In other words we are conservatively assuming no redshift evolution for 
the maximum value of k in the Fisher matrix analysis. This appears to be confirmed by A^-body simulations on the 
range0<z<2in|l74l|7a. 

It is not straightforward to derive a simple prescription for the expected uncertainties as a function of redshift and 
sky-coverage. Our simple analysis shows that high-redshift surveys (z > 1) covering a large fraction of the sky 
corresponding to a volume of about 100 Gpc^ might provide a 1-cr error on the running parameter of the order of 
A«NG - 0.4(50//i°L"). 



C. Constraints on equilateral non-Gaussianity from the galaxy bispectrum 



As shown by ||27| . the large correction to the galaxy power spectrum given by Eq. (IV.46i for the local model 
is essentially negligible for the equilateral one. On the other hand, higher-order correlation functions of the galaxy 



TABLE I: Characteristics of the galaxy surveys considered: sky coverage (Afl), redshift range {z„,i„ < z < z„ma). volume (V) 
and mean galaxy density (iig) together with the expected l-cr uncertainties on the amplitude and running parameters /ml and /inq 
from the Fisher matrix analysis of the galaxy power spectrum for the local model with fiducial values f^^ = 50 and «ng = and 
from the Fisher matrix analysis of the galaxy bispectrum (assuming the only the primordial non-Gaussian component, see text) /or 
the equilateral model with fiducial values = 100 and Mng = 0. Uncertainties are marginalized over the other non-Gaussian 
parameter, while cosmology and bias factors are assumed as known. 



Survey 


An [deg-] 




V[h-^Gpc^] 


iig [h^Mpc-^] 


A/ITl 






A«ng" 


Spectroscopic 


BOSS 


10, 000 


< z < 0.7 


5.8 


2.7 X 10^-* 


52 


2.6 


75 


1.5 


ADEPT 


28, 000 


1 <z<2 


114 


9.4 X 10-* 


9.7 


0.49 


22 


0.25 


Euclid 


20,000 


0<z<2 


108 


1.5 X 10^' 


84 


0.42 


17 


0.16 


HETDEX 


200 


2 < z < 4 


2.8 


5.0 X 10-* 


96 


3.8 


86 


0.60 


CIP 


300 


3.5 < z < 6.5 


3.6 


5.0 X 10-3 


61 


2.3 


35 


0.16 


Photometric 


LSST 


30, 000 


0.3 < z < 3.6 


390 


2.8 X 10-' 


3.5 


0.17 






PanSTARRS 


30,000 


< z < 1.2 


60.8 


1.7 X 10-3 


11 


0.59 







"Assumes a fiducial /^"f = 50. 
'Assumes a fiducial /j^J^ = 100. 
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FIG. 7: 1-cr contour plots for /nl and wng for the surveys considered for the local model from the power spectrum analysis. The left 
panel assumes an ideal all-sky survey from redshift z = to z = 3 and galaxy number density Hg - 10"^ Mpc"''. The right panel 
refers to the survey characteristics specified in Table |l] Notice that the choice of redshift bins of size Az = 0.5 and fixed minimal 
scale corresponding to = 0.03 /iMpc"' implies that the same range of scales is probed by all surveys, resulting in the same 
degree of degeneracy between the two non-Gaussian parameters /nl and hnq. 



distribution are expected to present an initial component when primordial pertiirbations are non-Gaussian. Previous 
analyses indicated that while current galaxy surveys are not yet competitive with CMB observations, in the future 
they could in principle provide an important confirmation of a possible CMB detection, 1 19] |20] [21]. These results, 
however, were based on the tree-level expression for the galaxy bispectrum in perturbation theory, which assumes 
the non-Gaussian contribution to the galaxy bispectrum to be given just by the primordial component to the matter 
bispectrum. In a couple of recent works, Jeong and Komatsu 1,37,1 and Sefusatti |38 1 have shown that non-linearities in 
the galaxy bias relation can induce significant 1-loop coiTections, corresponding to the leading contributions at large 
scales. In particular, such corrections depend on the effects of non-Gaussianity on the matter trispectrum. The effect 
on the galaxy, or halo, bispectrum is therefore similar to that on the power spectrum, where large-scale non-Gaussian 
corrections are due to the matter bispectrum. 

These preliminary results still require proper comparisons with A^-body simulations before being able to provide 
accurate predictions. For this reason, we make the very conservative choice to assume that the galaxy bispectrum is 
described by its tree-level expression in perturbation theory, where the only non-Gaussian effect is given by the initial 
matter bispectrum. In fact, it is reasonable to expect that the connections studied in |37, 38|-which we neglect — would 
significantly improve the constraints that we can place on primordial non-Gaussianity once reliable predictions for the 
bias parameters, derived for instance from numerical simulations, are given. 

We will therefore consider, for our Fisher matrix analysis, the reduced galaxy bispectrum Qg given by 



Qo(ki,k2,ki) = 



Bg(ki,k2,k2) 
Pg(ki)Pg(k2) + eye. 

^ [QcihMM) + QiihMM)] + % (IV.56) 



where Qc represents the reduced matter bispectrum induced by gravitational non-linearities, evaluated at second order 
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k [h Mpc"'] e 

FIG. 8: Effect of a running non-Gaussian component on the matter bispectrum. Left panel: reduced matter bispectrum for equilateral 
configurations, Q(k), as a function of the wavenumber k. Right panel: reduced matter bispectrum as a function of the angle 6 
between two fixed wavenumbers ki = 0.01 /jMpc"' and ^2 = 2fci. In both cases, the black continuous line indicates the Gaussian 
component while the red and blue continuous lines correspond to the full bispectrum including the non-Gaussian component for 
/nl = +100 and /ml - 100, respectively, at redshift z = 1. In addition, the dashed lines correspond to «ng = +0.5, dotted lines to 
"NG = -0.5 for both /nl = +100. 



in perturbation theory, Qi represents the reduced matter bispectrum due to non-Gaussian initial conditions and where 
the only effect of non-linear, quadratic bias is given by the constant term b2lb\. 

In Figure [8] we show the effect on the tree-level reduced matter bispectrum of equilateral non-Gaussian initial 
conditions. On the left panel the reduced bispectrum for equilateral configurations is plotted as a function of A; = = 
k2 - k^. The upper continuous curve coiTesponds to = 100 with «ng = 0, while the dashed and dotted curves 
corresponds to hng - 0.5 and hng = -0.5 respectively. The horizontal continuous curve represents the tree-level 
prediction for Gaussian initial conditions. The lower curves are for = -100 with «ng = and ±0.5. On the 
left panel we show the reduced bispectrum for configurations with sides ki - 0.01 /iMpc ' and ^2 - 0.02 /jMpc"' 
as a function of the angle between the vectors ki and k2. One can notice how the largest non-Gaussian corrections 
correspond to triangular configurations close to equilateral, while they vanish for squeezed configurations with 6 - Q 
and = n as expected for the equilateral model. 

For the details of the Fisher matrix analysis we refer the reader to Section III of BTll . since our calculation simply 
adds the extra running parameter and assumes a slightly different cosmology. We remind the reader, however, that 
unlike the case of the power spectrum, we consider here, for a given survey, all measurable triangular configurations 
with sides in the range defined by the fundamental frequency kf = In/V^^^ and a maximum wavenumber k„ia.x iden- 
tified by the non-linear scale which therefore increases with the median redshift of the survey or the specific survey 
redshift bin. This implies that different geometries and redshifts will probe different ranges in scales resulting in dif- 
ferent degeneracies between the two non-Gaussian parameters and n^c for our choice of the pivot point kp = 0.04 
Mpc"'. 

As for the power spectrum analysis of Section IV B we assume the bias parameters to be known and given by 
the halo model prediction plus a prescription for the halo occupation distribution (HOD) depending on the expected 
galaxy number density. While this might seem to be a strong assumption, we remark that an equally strong correlation 
between the linear and quadratic bias parameters b[ and b2 is expected in the halo model framework and observed 
in A^-body simulations and observations, see |88| and references therein. As shown in Sefusatti and Komatsu 121], 
if the halo model prediction for the halo bias functions, which can be measured in numerical simulations is assumed 
then the uncertainty on the non-Gaussian parameter marginalized over the HOD parameters is very close to the 
unmarginalized ones. A marginalization over the galaxy bias parameters bi and b2, in other words, would be an unduly 
conservative approach as it would introduce essentially unphysical degeneracies. 
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FIG. 9: 1-cr contour plots for and wng for the equilateral model from the galaxy bispectrum analysis. Left panel: ideal all-sky 
survey up to redshift z = 3 in redshift intervals of size Az = 1 and galaxy number density = 10"'/i^Mpc"\ Right panel: 
expectations for the surveys BOSS (outer dashed line), HETDEX {outer dotted line), CIP (continuous line), ADEPT (inner dotted 
line) and Euclid (inner dashed line). Note that the degeneracy between the two parameters depends on the choice of k„,„^ for each 
survey, as explained in the text. 



In Fig. [9] we show the results of the Fisher matrix analysis of the galaxy bispectrum in terms of the non-Gaussian 
parameters and hng- The left panel shows the l-cr contours for an ideal all-sky survey up to redshift z = 3 in 
redshift intervals of size Az - 1 and galaxy number density «g - 10"^ Mpc"^. Each interval assumes redshift bins 
of size Az = 0.5 with all quantities evaluated at the mean redshift of the bin and we ignore cross-correlations between 
different bins. In the right panels we instead show the results corresponding to the spectroscopic surveys of Table |l] 
where the marginalized constraints on and hng are given. 

It is interesting to notice how future redshift surveys such as Euclid will be able to place constraints on the running 
parameter of the order of Ahng - 0.3 at 95% C.L., if a significant non-Gaussian component of ^ 100 is detected. 
On the other hand, high-redshift surveys will be able to probe smaller scales and provide complementary information 
as shown by the different degeneracy in the /nl-wng plane for experiments such as HETDEX and CIP. We remind 
the reader that these simple results do not take into account the large-scale corrections to the galaxy bispectrum that 
come from the non-Gaussian effects on the matter trispectrum described in Il37ll38ll and that we expect to improve the 
constraints on both non-Gaussian parameters /nl and hng- A proper and complete assessment of the potential of LSS 
observations, including a joint analysis of power spectrum and bispectrum, requires a deeper and more quantitative 
understanding of effects which have just began to be studied and must necessarily wait for future work. 



V. CONSTRAINTS FROM JOINT CMB AND LSS OBSERVATIONS 



In this Section we consider the combination of the results from the CMB analysis of Section III with the results 
from the analysis of large-scale structure observables of the previous Section, both for the local and the equilateral 
model. In Table [ll] we show the expected 1-cr constraints on both non-Gaussian parameters /nl and hng for the local 
and equilateral model for all CMB and LSS experiments previously considered. We also compute the joint constraints 
between CMB and LSS for experiments providing results of the same order of magnitude. All errors are marginalized 
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over the other non-Gaussian parameters while we assume cosmological and bias parameters as known. All errors also 
assume the pivot point kp = 0.04 Mpc"'. In order to give a quantitative idea of how much the degeneracy between 
/nl and hng — which strongly depends on the choice of kp — is affecting the results, we include in parenthesis the 
corresponding unmarginalized errors. 

As a general remark, we expect large-volume and high-redshift galaxy surveys such as ADEPT, Euclid and LSST 
to provide constraints on the amplitude, as well as running of primordial non-Gaussianity comparable or even better 
than the constraints we expect from a CMB mission like Planck. In fact, from our simple analysis, we can expect a 
sensible improvement even over the results of the proposed CMBPol mission. 

In Fig. 10 we show the 1-cr uncertainties contours from the CMB bispectrum and LSS power spectrum analysis 
together with the combined constraints for the local model assuming the fiducial values = 50 and «ng = 0. In par- 
ticular, we compare the CMB experiments Planck and CMBPol with the LSS surveys Euclid and LSST together with 
the specific combinations of Planck with Euclid and LSST and CMBPol with LSST. The limits provided by Planck are 
comparable and complementary to those of LSST, the best large-scale structure probe in our choice. CMBPol would 
instead provide constraints quite close to the ideal case, only slightly improved by power spectrum measurements. 
We notice that although the choice of the pivot is not optimal for neither the CMB or the LSS observations as both 
individually show a significant degeneracy between the non-Gaussian parameters, in the combined constraint, such 
degeneracy is quite reduced. In particular, this shows the complementarity of the different probes when a running of 



TABLE II: Expected l-tr uncertainties on the amplitude and running parameters /nl and from the Fisher matrix analysis of the 
CMB bispectrum, the galaxy power spectrum for the local model with fiducial values = 50 and ;ing = and from the Fisher 
matrix analysis of the galaxy bispectrum (assuming only the primordial non-Gaussian component, see text) for the equilateral 
model with fiducial values f^ = 100 and Mng - 0, together with the corresponding joint constraints from CMB and large-scale 
structure observations. Uncertainties are marginalized over the other non-Gaussian parameter assuming the pivot point kp = 0.04 
Mpc"', while cosmology and bias factors are assumed as known. In parentheses we also provide the corresponding unmarginalized 
values. Joint constraints are shown only for the combinations which provide a sensible improvement over the CMB results. 
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CMB 



WMAP 


58.0 


(20.0) 


0.68 


(0.24) 


196 


(88) 


1.1 


(0.47) 


Planck 


4.7 


(3.1) 


0.10 


(0.07) 


30 


(29) 


0.30 


(0.29) 


CMBPol 


1.6 


(1.4) 


0.05 


(0.04) 


15 


(15) 


0.17 


(0.16) 


LSS, spectroscopic 


BOSS 


52 


(40) 


2.6 


(2.0) 


75 


(57) 


1.5 


(1.1) 


ADEPT 


9.7 


(8.7) 


0.49 


(0.43) 


22 


(12) 


0.25 


(0.14) 


Euclid 


8.4 


(7.3) 


0.42 


(0.37) 


16.7 


(6.7) 


0.16 


(0.06) 


HETDEX 


96 


(48) 


3.8 


(1.9) 


86 


(32) 


0.60 


(0.22) 


CIP 


61 


(30) 


2.3 


(1.1) 


35.2 


(8.9) 


0.16 


(0.04) 


LSS, photometric 


LSST 


3.5 


(3.1) 


0.17 


(0.15) 










PanSTARRS 


11 


(10) 


0.59 


(0.53) 
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(23) 


0.20 


(0.20) 


Planck + ADEPT 










16 


(15) 


0.17 


(0.16) 


Planck + Euclid 


3.6 


(2.8) 


0.09 


(0.07) 
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(7.3) 
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(0.07) 


Planck + LSST 
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(2.2) 


0.07 


(0.06) 










CMBPol + ADEPT 
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(11) 


0.13 


(0.12) 


CMBPol + Euclid 










10.6 


(6.7) 


0.10 


(0.07) 


CMBPol + Euclid + CIP 










9.5 


(5.7) 


0.07 


(0.04) 


CMBPol + LSST 


1.4 


(L3) 


0.04 


(0.04) 











"Assumes a fiducial /^"f = 50. 
'Assumes a fiducial /j^J^ = 100. 
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FIG. 10: 1-cr contour plots for /nl and n^Q for the local model for the Planck {dot-dashed blue) and CMBPol {continuous, blue) 
experiments compared to the LSS surveys Euclid {dotted, black) and LSST {dashed, black) together with the combined contours 
given by Planck and Euclid {thick, dotted, red), Planck and LSST {thick, dashed, red) and by CMBPol and LSST {thick, continuous, 
red). 



f^^' is considered. Notice that in the case of Planck, and even more so for CMBPol, such complementarity is a conse- 
quence of the different distribution of the non-Gaussian signal, for the specific local model, over the relevant range of 
scales, since the two observables are essentially probing the same range. Clearly the galaxy power spectrum, under the 
simple assumptions of our analysis provides a marginal improvement to the constraints from the ideal CMB results. It 
is reasonable to expect, however, that comparable or better constraints on local non-Gaussianity can be derived from 
an analysis of the galaxy bispectrum in light of the recent results of If37ll38l . which we did not consider here but that 
will be studied in future works. 



In Fig. 11 we show the 1-cr uncertainties contours from the CMB and LSS bispectrum analysis together with the 



combined constraints for the equilateral model assuming the fiducial values = 100 and hng = 0. The left panel 
shows the Planck experiments compared to the redshift surveys HETDEX, AI)EPT and Euclid, together with their 
combinations. The right panel shows instead the CMBPol against the surveys ADEPT, CIP and Euclid, together 
with their combined contours. In the case of the equilateral non-Gaussianity, future and proposed CMB and LSS 
experiments are expected to provide comparable constraints. In addition, they also tend to provide similar degeneracies 
between /nl and hng- 

This analysis essentially extends to the running parameter n^c the analysis of Sefusatti and Komatsu ETI on the con- 
straints on the amplitude of non-Gaussianities from the galaxy bispectrum in high-redshift surveys. Again, we should 
keep in mind that such analysis is quite conservative since we are neglecting large-scale non-Gaussian corrections due 
to non-linear bias 1.37., 38J . 



VI. DBI INFLATION: A CASE STUDY 



In this Section we present an example of how our analysis can translate into constraints of a specific model of 
inflation. As a case study we choose models in which the kinetic term is of the Dirac-Born-lnfeld (DBI) type instead 
of the canonical one. One of the virtues of such models is that they arise naturally in string theory |34| when inflation 
is driven by the motion of a D-brane li89J . Models of DBI inflation are particularly interesting in the present context 
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/nl /nl 

FIG. 11: l-cr contour plots for /nl and ;jng for the equilateral model from CMB and LSS measurements. Left panel: Planck (dot- 
dashed, blue) experiment compared to the LSS surveys HETDEX (dotted, black), ADEPT (dashed, black) and Euclid (continuous, 
black) together with the combined contours given by Planck with HETDEX (thick, dotted, red), with ADEPT (thick, dashed, red) 
and with Euclid (thick, continuous, red). Right panel: CMBPol (dot-dashed, blue) experiment compared to the LSS surveys CIP 
(dotted, black), ADEPT (dashed, black) and Euclid (continuous, black) together with the combined contours given by CMBPol 
with ADEPT (thick, dotted, red), with Euclid (thick, dashed, red) and with Euclid and CIP (thick, continuous, red). 



because they typically lead to large scale-dependent non-Gaussianity in the primordial spectrum of perturbations. We 
will first review the basics of DBI inflation and then possible concrete realizations. 

The dynamics of a D-brane in string theory is governed by the DBI action |90|. Consider a compactification from 
ten to four dimensions with a generic warp factor. For a spacetime-filling D3-brane the four-dimensional low-energy 
effective action is given by^ 



' DBI 



(VI.57) 



where /(0) is proportional to the warp factor at the position of the D3-brane in the compact space. The position is 
given by six scalar fields; in general this would lead to a multi-field inflationary model. In the explicit realizations we 
will consider, we assume that the brane moves along a single direction parameterized by the single scalar field 0. 

Homogeneous solutions for the above action coupled to gravity can be found 1341 [35 l| using the Hamilton-Jacobi 
formalism. One finds an upper bound for the velocity of the inflaton, (p < /(0) In the presence of strong waiping, 
f((f>) can become very large'' leading to a slow evolution that can support a prolonged stage of inflation in the presence 
of a slow-roll steep potential V(4>). In analogy with special relativity, it is convenient to introduce the function 

y=[l-m4>^l"\ (VI.58) 



^ Our convention for the metric signature is (-, +, +, +). 

^ f(<p) has the dimension of an inverse energy density. Here "large" means /m^ » 1, where m , is the string scale. 
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For flat potentials and small warping, y - 1, the dynamics of the system is the same as in slow-roll inflation. On 
the other hand, for steep potentials, i.e. a large slow-roll parameter rj > I, and strong warping we get y » 1. In this 



"relativistic" regime all higher-derivative terms in the expansion of the square root in Eq. VI.57 become important. 
The speed of sound of 0-perturbations becomes very small, c, = y"'. As first pointed out in |35| and calculated 
in detail in ll36l [58l . this can induce large non-Gaussianity in the primordial spectrum of density perturbations. The 
largest contribution comes from the equilateral configurations and can be estimated by^ 

/nl ^ -0.32 y2 = -0.32 cf . (VI.59) 

The running is straightforwardly calculated and is given by 

UNc s — ^ - -2^ = -2^ , (VI.60) 
a In K CsH 

where we defined s as the dimensionless time variation of the speed of sound^. 

The actual value of fNL and njvc depends on the model and in particular on the functions and f{(p) in Eq. VI.57 
These in turn depend on how the DBI inflation idea is specifically realized: in particular one has to specify how all the 
moduli are stabilized and the compact geometry in which the D3-brane moves. At present, no explicit realization of 
DBI inflation in string theory has been found that is consistent with the current CMB data. Rather than consider this 
to be a failure of the DBI inflation idea, we feel this reflects our limited knowledge of both moduli stabilization and 
the geometry of (compact) Calabi-Yau manifolds. 

In light of the above consideration, we deem it most sensible to consider effective models which are inspired by DBI 
inflation with phenomenological but realistic choices of /((/>) and V(^). The hope is that once a feasible embedding of 
DBI inflation in string theory is found, it might be reasonably approximated by the effective models that we analyze 
here. The string theory realizations proposed until now that have most closely reproduced existent data can be divided 
in three broad classes: UV models with > 0, IR models with /[(^(f)] < and angular models with f{<p{ty\ - 0. 

In the following, we describe each class and provide an archetypal choice of /(0) and V{(p). The resulting non-Gaussian 



parameters are plotted in Fig 12 together with the constraints from the data. 



A. The UV model 

In a broad class of models (initiated in ||9T1 ), the inflating D3-brane starts somewhere in the bulk of the compact 
dimensions and falls inside a region with increasingly stronger warping, referred to as the throat. As mentioned 
previously, we assume that the brane moves exclusively along the radial direction of the throat that we parameterize 
with (fi. We choose as zero for (p the tip of the throat where the warping is maximal. Inflation starts in the UV region 
(small warping) of the throat and ends somewhere close to the tip. Models in this class share the feature that /((/>) is 
monotonically increasing with time, as <p rolls towards zero. 

In the best-studied examples Il92ll93l . the throat is a conifold |94|. The resulting warp factor depends on how and 
where the singularity is cut off. A toy model that captures the features of a realistic solution {e.g. |92|) is given by 
ficj)) ^ A/((p^ + i?f, where A is a dimensionless parameter proportional to the D3-brane tension and yu determines the 
IR cut-off. Notice that the warp factor depends only on the radial and not on the angular position. Hence, a model of 



so-called angular inflation then takes place at constant /. We analyze this case in Section VI C 

As regards the potential, a series of detailed analyses ||9T1 195] |96l |97l |98l have confirmed that generically a mass 
term m^0^ is present, with m of the order of the Hubble scale. This fact constitutes a serious problem for slow-roll 
inflation (77 ~ (9(1)), but not for DBI inflation that proceeds with arbitrarily steep potentials. Other (/>-dependent terms 
can of course be present, but they will be subleading in the small-0 limit. 



' In I35II36| |58|. the sign convention of Maldacena |44| is used which is opposite to the WMAP convention. We adopt the latter, which explains 

the minus sign in Eq. |VI.59| 
* This is in analogy with the slow-roll parameter r/ ^ i/ {^H) 
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In 1991 the simplest realistic UV model of DBI inflation was compared with the data. It was found that no model 
exists that obeys all the constraints imposed by the internal consistency of string theory and that is compatible with 
the data. In IIIOOL the analysis was generalized to include other plausible f{(p) and V((p), but resulted in the same 
conclusion. As pointed out in IIIOII . the reason for these negative results can be drawn back to the existence of a 
geometrical limit for the allowed range of (/>. The upper bound is given by (Acp/Mpi)^ < AjN, where is a quantum 
number that characterizes the strength of the warping and for consistency » 1. 

If one could relax this bound the cosmological data could be easily reproduced [W, 100, '102 |. Several ideas and 
speculations about how this could be achieved in string theory have been proposed. Nevertheless, until now no explicit 
model has emerged that is fully under control and where the geometric bound is violated. 

Given all the above considerations, we focus on a string theory-inspired model of DBI inflation, for which we 
neglect the geometrical bound and select the simplest realistic functions /((^) and y(0). Our archetypal UV DBI 
model is chosen to be defined by 

/(-^^ = t:^;^^-^ ' (VI.61) 

which is known to induce a power-law period of inflation f34l. When (p 's> fi the warp factor approximates that of AdS, 
f((p) ~ A<p^'^. Defining ^ m^A/(6M^i) and assuming » 1, we find that the scale factor evolves as a{t) oc t'\ In 
ll36l . the /nl parameter produced with the above choices was estimated for the equilateral configurations. The result 
in the WMAP notation is 

(VI.62) 

where (pt indicates the value of around which the CMB perturbations exit the horizon. The running of the non- 
Gaussianity riffc is straightforwardly calculated as^ 

_ gln|/NLl 4 

"NG = ^, , - - ■ (VI.63) 

A few comments are in order. First, uj^q > 0, which implies that the speed of sound decreases as time evolves. This 
is because in the UV models the brane advances towards regions with larger warping, i.e. f > 0. Second, we notice 
that for simplicity we have chosen an archetypal model with three parameters '° A, /u and m. In the present context, one 
parameter is fixed by the COBE normalization while the other two still allow for enough freedom to cover most of the 



interesting region in the sector /^^ < n rii^Q > 0. For illustrative purposes, we plot in Fig. 12 the prediction of the 
model in the /nl-wng plane for different (constant) values of (pt/Mpi together with the expected 1-cr constraints from 
CMB experiments (left panel) and combined CMB and large-scale structure constraints (right panel). This clearly 
shows that the geometric bound on the range of (f> IIOII must be evaded in order to reconcile UV DBI models with 
experiments. 



B. The IR model 



In a second class of models (initiated in |103|), the inflating D3-brane starts at the tip of a throat (the IR region), 
where the warp factor is maximal, and evolves towards the bulk of the compact space. This motion still takes place 
along the radial direction, but opposite to the UV models since /(0) is now monotonically decreasing with time. 



11.12 



with respect to I36II58I where a scale invariant /nl would correspond to jing = 1- 



' We use a different convention (see Eq. ] 
From an effective field theory point of view, there would be another additional parameter, i.e. the reheating temperature T^, that determines the 
number of e-foldings between the end of inflation and the time at which the CMB perturbations have exited the horizon. This parameter would 
be degenerate with say /i in the present model. On the other hand, once an explicit embedding of both inflation and the standard model in a string 
theory construction is specified, Tr should be determined. 
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FIG. 12: 1-cr contour plots for CMB alone (left) and combined CMB and LSS measurements (right) vs. DBl predictions. Assumes 
fiducial values f^j^ = -50 and = 0, the first corresponding to the lower 1-cr limit from the WMAP bispectrum analysis 1321 . 



In flQA] an extensive comparison of this class of models with the data has been performed. Similar to the UV 
models, no IR DBI model was found which is compatible with the data and which obeys all the consistency conditions 
imposed by the embedding into string theory. The main obstacle is that the largest allowed number of e-foldings 
during the DBI regime is smaller than order ten. Even if inflation can be prolonged by a phase of fast or slow roll, the 
perturbations produced during the relativistic DBI regime are too red-tilted. References Il36lll04lll05l concluded that 
the model might be saved by assuming that the CMB perturbations are produced in a phase preceding the DBI ones. 
This would be a Hagedorn phase (i.e. the Hubble temperature is larger than the warped string scale) in which the string 
partition function formally diverges and all string states are infinitely populated. This scenario is too poorly understood 
to allow for a rigorous treatment of the production of perturbations, and so we prefer to neglect this possibility until a 
better description is available. 

As we did for the UV models, we consider an archetypal effective model of IR DBI inflation. To get a model 
compatible with current data, we neglect some of the constraints present in the known string theory constructions. In 
particular we relax the bound on the inflaton field range IIIOII and allow for an arbitrary number of e-foldings in the 
ultra-relativistic DBI regime. We choose 

/(•^^ = T:^;^^-^ ' V(cf>)^Vo-lmY, (VI.64) 

and we will often use the reparameterization = /3H^. As we mentioned above, considering various corrections 
[91 95j .96j .97. 98J one finds generically /3 ~ 0(1). In li36J . the /nl parameter produced by the above model in the 
equilateral configurations was estimated to be 

where A^, is the number of e-foldings prior the end of inflation around which the CMB perturbations are produced. 
The running of the non-Gaussianity was estimated to be 



(VI.66) 
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In contrast to the UV model, n^c < 0, a consequence of the fact that / < 0. Notice that for a large class of potentials, 
the sign of n^c does not depend on V{4>)\ as we will soon see, this is not the case in the angular model. In addition, our 
archetypal effective model has enough parameters to cover most of the interesting region in the sector f^i < n umg < 



0. We show in Fig. 12 some rif^cifNi) curves for various values of The curves are similar to those in the UV model 



but with opposite sign. 



C. The angular model 



We now consider a third model where f{(p) is constant, which in a sense lies in between the UV and IR models. In 
a string theory construction, this can be achieved if the inflationary D3-brane moves along the angular directions of 
a warped conifold, and so maintains a constant radius. An explicit realization of this was constructed in [1061 (see 
also 1 107 1 for an earlier attempt). The D3-brane moves along the tip (which is an 5^) of a warped deformed conifold, 
hence the name inflation at the tip. The motion is driven by the F-term potential generated by the presence of a stack 
of D7-branes with gaugino condensation. In this Section the field refers to the angular position of the D3-brane, 
while the radial position is constant. 

The analysis of 1 106 1 shows that, even for a mildly warped throat, a phase of ultra-relativistic DBI inflation generally 
occurs. On the other hand, once the bounds induced by the consistency of the string theory construction are imposed, 
not enough e-foldings of DBI inflation can be generated". Although some ideas have been proposed to overcome 
these difficulties 11061 . no successful concrete model of angular DBI inflation has been found. 

Analogously to the IR and UV cases, in the lack of a genuine string theory construction we consider an archetypal 
elTective model. For simplicity, we propose the most minimal model with just two parameters,'^ 

f{<P) = const , y(0) = ^m^^ . (VI.67) 
The resulting non-Gaussianity can be estimated using Eq. |VI.59] and the approximated solution obtained in Appendix 

El 

/NL--0.3/~-^mV. (VI.68) 
The running of /nl has been calculated in Appendix [C] The result is 

«NG--— (VI.69) 

where again (p - (p^ corresponds to the time around which CMB perturbations exit the horizon. If we use COBE 
normalization to fix one of the two parameters / or m and express the other in terms of /nl, we find that «ng(/nl) ~ 
10""^, i.e. it is constant at leading order. A small «ng could have been anticipated: hng takes the opposite sign in the 
UV and IR models, respectively, and the angular model lies in the middle, i.e. fuv > fang = > fjR, close to where 
njVG passes through zero. As can be seen from Fig. [12] such a small n^c is much below the experimental sensitivity. 
Contrary to the UV and IR models, in the angular model the sign of nj^c is not fixed by f{(p) but it depends on the 
potential. In particular, a convex V{(p) gives a negative hng and vice versa. 

VII. CONCLUSIONS 

We studied the possibility of constraining a scale-dependence of the amplitude of a non-Gaussian component in the 
primordial perturbations in future observations. Focusing on two models — local and equilateral — for the bispectrum 



' ' On the contrary, if one allows for fine tuning, models of slow-roll inflation can be found that are compatible with all present data. 

In this model N, ~ y(ip,/2)- which shows that a superplanckian field range is required to get enough e-foldings. Hence, a string theory 
construction needs to go beyond this minimal model. A problematic field range bound is therefore common to all three models of DBI inflation. 
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of the curvature fluctuations at early times, we performed a simple Fisher matrix analysis of the CMB bispectrum in 
terms of the non-Gaussian amplitude parameter /nl and a running parameter hng- In addition, we considered a simple 
Fisher matrix analysis of the large-scale structure power spectrum for the local model and of the galaxy bispectrum 
for the equilateral model of non-Gaussianity. We compared and combined these results to those expected from CMB 
observation to explore possible complementarities between such different probes of the early Universe. 

We find that, in the event of a detection of a relative large non-Gaussian component, upcoming and future CMB 
missions such as Planck and CMBPol are capable of providing significant constraints on a possible scale-dependence 
such as predicted in certain inflationary models like DBl inflation. Assuming the current central values of the limits 
from WMAP observations for the local and equilateral models, that is f^^- = 38 Q and =51 ||32l , our Fisher 
matrix analysis indicates that we could expect 1-cr uncertainties on the running parameters of the order of Ahng - 0. 15 
and A«NG - 0.6 from the Planck mission. The proposed CMBPol mission should reach uncertainties of the order of 
Amng - 0.05 and A«ng - 0.35 for the local and equilateral model respectively, quite close to those expected from the 
ideal CMB mission. 

We remark that the constraint on the extra running parameter come, in principle, at no cost with respect to the 
constraints on the amplitude. In fact, given a specific observable and a specific model it is always possible to choose 
the optimal pivot point that removes any degeneracy between the two parameters. However, it is useful to determine an 
optimal pivot point that can be used by different probes. We notice that, assuming our definition for the scale-dependent 
non-Gaussian parameter fj^hik) given in Eq. (11.7 1, for an experiment like the proposed CMBPol mission, quite close 
to the ideal CMB experiment, the optimal pivot point is about kp ^ 0.3 Mpc ' for the local model and kp ^ 0.4 Mpc ' 
for the equilateral model, indicating the different distribution of the signal across the range of observable scales of the 
two models. 

We also studied the constraints achievable by future large-scale structure surveys, finding results comparable to 
those from CMB measurements both for local and equilateral non-Gaussianity. In particular, for the local model we 
focused on the effect of non-Gaussian initial conditions on the bias of halos and galaxies which received considerable 
attention in the recent literature, 1 22 , 23 , 24 , 25 , 26 , 27 1 and which lead to constraints on the local f^^- from current 
observations already comparable to those from the CMB, Ii24l l26l . We considered an expression for the correction 
to the galaxy linear bias parameter given in terms of the initial matter bispectrum, as those proposed in ||231 l27l . 
which naturally allows for a scale-dependent /nl(^i , ^3)^ and derived the expected constraints on the non-Gaussian 
parameters /nl and hng from measurements of the galaxy power spectrum in future spectroscopic and photometric 
surveys. We find that large-volume surveys such as LSST and Euclid should be able to provide constraints smaller than 
those expected from Planck and somehow able to improve those expected from the proposed CMBPol mission both 
on /nl as well as on the running parameter «ng^ possibly achieving a 1-cr uncertainty below Ahng - 0.05 in a joint 
analysis. Moreover, we noticed that the CMB and LSS constraints show a marked complementarity in the specific 
case of local non-Gaussianity. 

The effects of equilateral non-Gaussianity on the galaxy power spectrum, on the other hand, are expected on theo- 
retical grounds, to be negligible, 1271 . In this case, we turned our attention to the galaxy bispectrum. Recent studies 
of the effects of non-Gaussian initial conditions on the bispectrum of biased populations have shown that large-scale 
corrections, similar to the one affecting the power spectrum in the local case, are present both for local and equilateral 
non-Gaussianity ll37l[38l . However, these preliminary results still need to be properly tested against numerical simu- 
lations before providing reliable predictions for the galaxy bispectrum. We therefore considered the simple tree-level 
expression for the galaxy bispectrum in perturbation theory where the only non-Gaussian correction is represented by 
the primordial component to the matter bispectrum. This results in a quite conservative Fisher matrix analysis since the 
neglected (1-loop) corrections due to non-linear bias are expected to significantly increase the non-Gaussian signature 
on the galaxy bispectrum, in principle for any non-Gaussian model. Nevertheless, our analysis, which constitutes an 
extension to the case of running non-Gaussianity of the analysis of 11211 . shows that future large spectroscopic surveys 
like Euclid could provide constraints on the hng parameter of the order of A«ng - 0.1, comparable or even better than 
those from CMBPol, while presenting a sensible complementarity. 

It is to be stressed that these simple and preliminary results do not constitutes a proper analysis of the ability of large- 
scale structure observations to constrain non-Gaussianity. They show, however, that the large-scale galaxy distribution 
can be an as sensitive, although less direct, probe of the early Universe and that a possible detection of a significant 
level of non-Gaussianity in the CMB should find a confirmation in the large-scale structure. Moreover, they indicate 
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that extra information, as a possible scale-dependence, on such non-Gaussian component can be extracted by a joint 
analysis of CMB and LSS data. Clearly a complete analysis of large-scale galaxy correlators should take into account 
the power spectrum as well as the bispectrum and their covariance, which we neglected here. Still, we expect our 
results to motive further studies in the same direction. 

Finally we considered a specific, string-motivated, inflationary model — DBI inflation — as a case study for running 
non-Gaussianity. This model naturally predicts a large non-Gaussian signature with a significant scale-dependence. 
Interestingly, for two common DBI inflaton Lagrangians, a smaller non-Gaussian amplitude corresponds to a stronger 
running. In this respect, constraints on this specific model should consistently take into account the peculiar scale- 
dependence of the predicted curvature bispectrum. We showed, in fact, that the determination of the running parameter 
fiNG can place additional constraints on the parameters of the specific model of inflation breaking degeneracies other- 
wise present from a simple determination of the overall amplitude of the primordial bispectrum. Although an explicit 
realization of DBI inflation in string theory has not been found, under a purely phenomenological point of view, 
our analysis shows how the extra information provided by the running parameter n^Q directly translates into extra 
information on the parameters of the model. 
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APPENDIX A: SCALE-DEPENDENCE AND TRIANGULAR CONFIGURATIONS 



In this Appendix we compare the results from the all-sky Fisher matrix analyses of the galaxy power spectrum 
(for local non-Gaussianity) and of the galaxy bispectrum (for equilateral non-Gaussianity) obtained assuming the 
geometric mean, adopted in this work to define the scale-dependent non-Gaussian parameter /nl(^i , fe, = /nl(^), 
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FIG. 13: Correlation coefficient ci2 between the two non-Gaussian parameters /ml and «ng as a function of /„,„^ assuming several 
values of the pivot point kp in the T + E case for CMBPol. Left panel: local model. Right panel: equilateral model. 



given by 

with the same results obtained instead assuming the arithmetic mean 



(A.l) 



(A.2) 



assumed by other works and derived as the "coiTect" form of the curvature bispectrum predicted by DBI inflation. In 
Fig. 13 the thick lines coiTespond to the geometric choice of Eq. while the thin lines to the arithmetic one of 
Eq. (A.2 1. On the left panel, we show the analysis for the galaxy power spectrum and local non-Gaussianity, while 
on the right panel we show the analysis for the galaxy bispectrum and equilateral non-Gaussianity. In the case of 
equilateral non-Gaussianity, the difference is minimal (and, as expected, almost unnoticeable) since most of the signal 
is coming from equilateral configurations. The two choices show instead a significant diff'erence in the case of local 
non-Gaussianity, both in terms of degeneracy and overall marginalized errors. While the factorizable, geometric mean 
of Eq. (A.l I is motivated by the easier implementation of the CMB estimators, we can assume this choice to have a 
purely phenomenological value. On the other hand, in the equilateral model, we can expect the CMB results to hold 
even assuming the arithmetic mean, better motivated under the theoretical point-of-view. A proper implementation 
of an estimator for the CMB bispectrum described by a running non-Gaussian parameter fniXK) with K given by 
Eq. (A.2 1 will be considered in a future work. 



APPENDIX B: THE CHOICE OF THE PIVOT POINT 



The degree of degeneracy between the non-Gaussian parameters /nl and hng strongly depends on the specific 
choice for the pivot point the enters the definition of Eq. pi.7| i. Since such a choice is arbitrary, it is convenient to 
choose a value of kp that minimizes the correlation between the parameters. However, the best value obtained in this 
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FIG. 14: Correlation coefficient cn between the two non-Gaussian parameters /nl and /ing (see text) as a function of the pivot point 
kp for the temperature (thin lines) and temperature and polarization analysis (thick lines) assuming the Planck (dashed lines) and 
CMBPol (continuous lines) experiments with l„,ax = 2500. Left panel: local model. Right panel: equilateral model. 



way depends in turn on the observable considered and on its characteristic signal-to-noise over the range of scales 
probed, which, in addition, is specific to the non-Gaussian model under consideration. Since the CMB bispectrum is 
currently providing the best constraints on primordial non-Gaussianity and also constitutes its best-studied probe to 
date, in this work we determine the value of pivot point kp to be the value that minimizes the degeneracy between /nl 
and hng as derived from the CMB analysis. 

Given the CMB Fisher matrix Fjj, the correlation coefficient for the parameters is given by 

cn^ , ^ (B.l) 

where the indices 1 and 2 label the parameters /nl and hng. respectively. Notice that since (F"')i2 ~ f^^l^, (/^"')22 ~ 
/j^^ ^ii'l ')ii does not depend on /nl, the correlation coefficient is independent on the fiducial value of /nl chosen. 

In Fig.[T4]we plot the correlation coefficient between the two non-Gaussian parameters /nl and hng as a function 
of the pivot point kp for the Planck experiment (dashed, blue lines) and the CMBPol experiment (continuous, red lines) 
assuming in both cases v - 2500. The thin lines correspond to the temperature analysis while thick lines correspond 
to the combined temperature and polarization information. Finally, on the left panel we assume local non-Gaussianity 
while on the right panel equilateral non-Gaussianity. As can be expected, the local model prefers a lower value for 
kp of about 0.03 Mpc ' where ^0.1 for CMBPol, due to the higher signal-to-noise in configurations correlating 
small and large angular scales. For the equilateral model we find instead a best value around kp - 0.05 Mpc where 
c\2 - 0.06. The corresponding values for Planck are naturally lower in both cases. 

A high value of kp can be expected to be a better choice for an analysis of large-scale structure information, since 
the galaxy distribution and even more cluster abundance probe smaller scales. For simplicity we could choose to 
determine kp as the optimal value for the CMB bispectrum analysis, so that the expected errors on /nl from future CMB 
experiments, marginalized over the running parameter n^c, are not significantly diff'erent from the values already 
determined in the literature under the assumption of a constant /jl. In this respect the choice of kp ^ 0.04 Mpc"', 
proposed already by Lo Verde et al. |13|, seems to be a particularly reasonable choice as a compromise between the 
two models of non-Gaussianity and in the perspective of future CMB missions beyond Planck. 
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APPENDIX C: RUNNING OF NON-GAUSSIANITY IN ANGULAR DBI INFLATION 



In this Appendix we obtain an analytical estimate for the running of non-Gaussianity in the simplest model of 
angular DBI inflation, i.e. with constant / and V - m^(p^ 12. Our starting point is 

HiVc - -2i = -2-^ = 2^ ^ 2 — ^— -, (C.l) 
c,H yH H'H^U 

where in the last step we used the fact that we are in the ultra-relativistic regime 

r = ^l+4/(//')2 - 2//' V7 » 1 ■ (C.2) 
Using the Hamilton- Jacobi formalism it is straightforward to solve for H at linear order in 

md) m 

H{cf>) ^ ^ = —{cf,. - Acf>) , (C.3) 

where we have introduced the distance A(f> traveled by the inflaton since the CMB perturbations have exited the horizon 
at = 0t. To obtain the running we need to go beyond the linear approximation. We therefore make the ansatz 

H{(P) = + ho] + (-^ + hi] A(f) + h2A(f>^ . (C.4) 

We now expand 



V6 / I V6 



3//2 = y + (C.5) 



at second order in A0, using the exact expression for y given in Eq. C.2 Then we linearize in /zq, hi and /z2 and solve 
for them. The result is 

2 1 

4>^hQ = (pthi = /l2 = ;=-r . (C.6) 

Plugging this solution into Eq. |C.l| we get the estimate for the running of non-Gaussianity produced close to (f> = <pt 

8 

nNG ^ - o ■ (C-7) 
We have numerically checked that this is a good estimate for y » 1 and A(p I. 



